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Topological entanglement entropy, a measure of the long-ranged entanglement, is related to 
the degeneracy of the ground state on a higher genus surface. The exact relation depends on 
the details of the topological theory. We consider a class of holographic models where such 
relation might be similar to the one exhibited by Chern-Simons theory in a certain large 
N limit. Both the non-vanishing topological entanglement entropy and the ground state 
degeneracy in these holographic models are consequences of the topological Gauss-Bonnet 
term in the dual gravitational description. A soft wall holographic model of conhnement 
is used to generate hnite correlation length but keep the disk topology of the entangling 
surface in the bulk, necessary for nonvanishing topological entanglement entropy. 
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1 . Introduction and summary 

Entanglement entropy has recently emerged as an important qnantity whose signih- 
cance spans varions snbjects ranging from qnantnm gravity to qnantnm Hall effect. One 
of the properties of the latter is nonvanishing topological entanglement entropy, which can 
be defined to be a finite term 7 in a large radins expansion of the entanglement entropy 
for a disk region in a 2+1 dimensional theory with finite correlation length [ 1 , 2 ]: 

See = oiR — 1 + ■ ■ ■ ( 1 - 1 ) 

More generally, topological entanglement entropy can be a good non-local order parameter 
for qnantnm liqnids with long-range order in sitnations where more conventional local order 
parameters are not nsefnl. 

Many qnantnm Hall systems can be described by Chern-Simons theories, where topo¬ 
logical entanglement entropy can be compnted; it eqnals the Sq component of the modnlar 
S-matrix of the theory [3]. It can also be shown to represent a constant term in a par¬ 
tition fnnction on a sphere or a logarithm of the total qnantnm dimension. Topological 
entanglement entropy is also related to the degeneracy of the gronnd states for theories 
compactified on a spatial snrface of genns g. In particnlar, for abelian Chern-Simons 
theories, the relation is 

2177 = S, , ( 1 . 2 ) 

where Sg is the logarithm of the nnmber of the gronnd states. For a non-abelian Chern- 
Simons theory, the relation between Sg and 7 can be more complicated. As we review 
below, an interesting simplification occnrs for the SU{N)k Chern-Simons theory in the 
limit A7 3 > /c 3 > 1 . In this case, the relation between Sg and 7 is 

2(9-1)7 = 55. ( 1 , 3 ) 

Recently, a seminal work by Ryn and Takayanagi [4] has inspired work on entan¬ 
glement entropy in the context of AdS/CFT correspondence. Althongh qnite difficnlt to 
calcnlate by conventional field theoretic techniqnes, the entanglement entropy can be easily 
obtained in the dnal holographic description by compnting the minimal area of the bnlk 
snrface anchored on the entangling snrface at the AdS bonndary. This raises a natnral 
qnestion whether one can constrnct holographic theories with nonvanishing topological 
entanglement entropy 7 . Compnting 7 involves finding an extremal snrface in the bnlk 
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anchored to a large circle at the boundary of asymptotically AdS space. In [5] this excersise 
was performed for the AdS-soliton geometry, which describes a conhning 2+1 dimensional 
held theory. In that model, at large radii, the dominant bulk hypersurface has the topology 
of a cylinder, and gives rise to the vanishing topological entanglement entropy in accord 
with our expectations: such a QCD-like theory is not expected to have a ground state 
with nontrivial long-range order. For other work on topological entanglement entropy in 
the context of holography, see [6,7]. 

In this paper we show that non-vanishing topological entanglement entropy arises 
naturally for certain held theories whose holographic dual four-dimensional descriptions 
include a Gauss-Bonnet term. In four-dimensional gravity, the Gauss-Bonnet term is purely 
topological. Its contribution to the holographic entanglement entropy is also topological 
and can be computed by integrating the Euler density over the minimal surface in the 
bulk [8-10]. It produces the Euler characteristic of that surface and hence, to get a non¬ 
vanishing 7 , we need the minimal surface to have the disk topology in the bulk, rather than 
the cylinder topology. This is precisely what did not happen in the example studied in [5]. 
There are geometries where only a disk-like surface is a solution (and a cylinder-like is not) 
- see e.g. [11]. The model considered in [11] does not quite work for us, since it does not 
have a gap in the excitation spectrum, but it shows that we may have a chance by making 
the co nfi nement “softer”. Indeed, in this paper we show that certain soft-wall models 
of co nfi nement do support entangling surfaces with the disk topology and, moreover, a 
constant term in the entanglement entropy due to the Einstein-Hilbert part of the bulk 
action is absent. Hence, the addition of the Gauss-Bonnet term to these models ensures 
nonvanishing 7 . 

These four-dimensional holographic models (reviewed below) satisfy two important 
conditions: 

(i) The extremal surface anchored on a large circle at the AdS boundary has the topology 

of a disk. 

(ii) The constant term in the large radius expansion of the area of this bulk surface 

vanishes. 

The second condition comes from the expectation for 7 to come entirely from the Gauss- 
Bonnet term in the action, in the anticipation of its relation to the degeneracy of states 
on a genus g surface. This is because contribution of the Gauss-Bonnet term to Sg is also 
topological, while a contribution from the Einstein-Hilber term to the entropy is necessarily 
extensive (equals the area of the horizon). Indeed, we hud that for the soft-wall models 
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the relation between 7 and Sg is the same as for the SU{N)k Chern-Simons theories in 
the 3> /c 3> 1 limit, (1.3). To be more precise, this relation is only trne for the part of 
the gronnd state entropy which comes from the Ganss-Bonnet term. As it tnrns ont, in 
Einstein-Hilbert gravity it is very hard to End a holographic geometry, whose bonndary is 
a genns g > 1 snrface, with vanishing entropy at zero temperatnre. This statement is trne 
for pnre AdS and remains trne for models we consider. This finite area of the horizon term 
spoils the relation (1.3). Snch a term however is exponentially snppressed as the prodnct 
of the conhning scale and the size of the genns g snrface becomes large. 

The rest of the paper is organised as follows. In Section 2 we review some facts 
abont Chern-Simons theories. In particnlar we review relations (1.2) and (1.3). In Section 
3 we review some basic facts abont holographic fonr-dimensional Ganss-Bonnet gravity 
and point ont that a pntative held theoretic dnal of this gravitational model may possess 
non vanishing topological entanglement entropy and gronnd state degeneracy, related by 
(1.3). In Section 4 we constrnct holographic models which have non-vanishing topological 
entanglement entropy - the soft-wall models discnssed above. In Section 5 we stndy the 
gronnd state entropy of these soft-wall models compactihed on Riemann snrfaces of genns 
g. We discnss onr resnlts in Section 6 . The appendices A and B discnss the spectrnm 
and the entanglement entropy of a slab region in the soft-wall model. We also discnss the 
topological entanglement entropy in a certain holographic model of qnantnm Hall effect in 
appendix C. 


2 . Quantum dimensions and ground state degeneracy in Chern-Simons theory 

In this section we briefly review the physical interpretation of the qnantnm dimension 
appearing in the Chern-Simons theory and its relation to the entanglement entropy. We 
hnish the section by ontlining how to obtain the relation between topological entanglement 
entropy and gronnd state degeneracy on Tig x . 

It is a well known fact that Chern-Simons theory is related to the corresponding Wess- 
Znmino-Witten (WZW) theory [12]^. Consider the theory on T^. Under the modnlar 
S transformation, S : r “I/'?', exchanging two non-contractable cycles, the WZW 
character is transformed as 


X„(-l/r) = ^5^ Xb(r) . (2.1) 

b 


^ For an extensive review of WZW theory see e.g. [13]. 
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The modular S-matrix, plays a crucial role iu the correspoudiug Cheru-Simous theory 
siuce the fusiou rules [14] (see also [15] for a review iu coudeused matter applicatiou) aud 
partitiou fuuctious ou various Riemauu surfaces [12] cau be writteu as combiuatious of 
these matrices. 

We shall focus ou the quautity called quantum dimension^ da = Sa/Sg. This object 
has a special iuterpretatiou as a relative degeueracy of state iu represeutatiou a compared 
to the degeueracy of state iu the ideutity represeutatiou, deuoted as a = 0 [16]. This cau be 
showu more precisely by realisiug that the degeueracy of represeutatiou a cau be writteu 
Cl'S Xa(q 1), where q = Thus, we have 


da = lim 


Xa(q) 


= lim 


Eb-^aXbig) 


* 50 ° 


The topological eutauglemeut eutropy iu (1.1) cau be writteu as (se e.g. [3]) 


7 = - log 5'o = log V 


V = 



( 2 . 2 ) 


(2.3) 


where V is called the total quantum dimensions [1,2]. Recall also that Sg is the partitiou 
fuuctiou of the Cheru-Simous theory ou [ 12 ]. 

We proceed by poiutiug out the relatiou betweeu 7 aud Sg where Sg is the eutropy of 
the Cheru-Simous theory ou Eg X 

Sg=\OgZ[EgXS^], (2.4) 

where a simple relatiou ( 1 . 2 ) cau be fouud for U{l)k gauge group aud (1.3) for SU{N)k 
with iV > /c > 1 

2.1. U{l)k abelian Chern-Simons theory 

For abeliau Cheru-Simous theories, the elemeuts of the S-matrix are just phases aud, 
cousequeutly, the total quautum dimeusiou is simply 


V = 


^ |da|2 = Ek , 


\ a=l 


(2.5) 


( 2 . 6 ) 


aud, therefore, the topological eutauglemeut eutropy for this theory is 

7 = log T> = 1 log(/c) , 

There are mauy ways to get the degeueracy of the 17(1)fc theory ou a geuus g surface. 
For the least mathematically iuvolved approach, see [17]. The grouud state degeueracy 
aud the associated eutropy are 

Z[ZgXS^-U{l)k] = k3 ; Sg = g\ogk. (2.7) 

Thus, the topological eutauglemeut eutropy aud grouud state eutropy are related by (1.2). 
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2.2. SU{N)k Chern-Simons theory 

We now consider the non-abelian Chern-Simons theory with the SU{N)k gange gronp 
where 3> /c. The expression for the total qnantnm dimension can be fonnd in e.g. [3] 
bnt the calcnlation for large N limit can be qnite involved. For an alternative method 
of compnting total qnantnm dimensions, see e.g. [18,19]. A simple way to hnd a nice 
expression for the limit we are interested is to nse the level-rank dnality i.e. the partition 
fnnction of SU{N)k and SU{k)N on are related by [ 20 ] 

Z[S^ SU(N),] pk 

Z[S^,SU{k)t,] VN' ’ 

One can nse the expression for Z[S^, SU{k)N] for hnite k and N ^ oo, which is also 
presented in [ 20 ], 

log Z[S^,SUik)N] ^ —{k^ -1) log N + 0{N^) . (2.9) 

Using (2.8) and (2.9), the total qnantnm dimension T> in this limit is 

2 , A’2 

p^jY+fc /2 . 7 = logT> = — logiV . (2.10) 

As for the gronnd state degeneracy on Eg x we nse the same approach ontlined 
above. The level-rank dnality for this manifold is fonnd in [21] to be 

Z[Z,xS\SU{N),] fNy 

Z[ZgX S\SU{k)N] \k J ' ’ 

The expression for Z[Zg X S\SU{k)N] with N ^ k can be fonnd in [22] 

logZ[Zg X S\ SUik)N] l)ik^ - 1) logiV + 0{N^) . (2.12) 

As a resnlt, the entropy on Eg x with the gange gronp SU{N)k, can be expressed as 
follows: 

Sg = logZ[Zg X S\SU{Nk)] =glogiN/k) + (g - l)ie - 1) logiV + O(iVO) . (2.13) 

In the limit 3> /c 3> 1 , we can then relate 7 and Sg for SU{N)k nsing ( 2 . 10 ) and (2.13), 
and the relation between them becomes (1.3). 
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3. Gauss-Bonnet holography in 4 dimensions, entanglement entropy and 
ground state degeneracy on Eg x 

The Gauss-Bonnet theory is one of the simplest extensions of the Einstein gravity. It 
is described by the Einstein-Hilbert action with a 4-dimensional Euler density, the Gauss- 
Bonnet term, added. The action of this theory is 


where 








(3,1) 


(3.2) 


The Gauss-Bonnet term is non-dynamical, but affects physical quantities, such as the en¬ 
tropy and the entanglement entropy. The correction term is simply an additional constant 
proportional to A. The Gauss-Bonnet gravity is problematic in the following ways. It has 
been shown that for the positive Gauss-Bonnet coupling, A > 0, one can merge black holes 
and violate the second law of thermodynamics [23]. We will mostly be interested in the 
situation with the negative coupling, A < 0, where one can have black holes with negative 
entropy [24,25,26,27]. Also, graviton scattering in higher-dimensional Gauss-Bonnet the¬ 
ories exhibits violation of causality [28]. Nevertheless, none of these issues will appear in 
the present work. 

In this section, we review computations of the entanglement entropy and the black 
hole entropy in the presence of the Gauss-Bonnet coupling. We show that contributions 
from the Gauss-Bonnet term to the ground state entropy and to the entanglement entropy 
are related by (1.3). 


3.1. Entanglement entropy and the entropy of topological black hole 


To calculate entanglement entropy for the held theory dual to 4-dimensional Gauss- 
Bonnet gravity, one has to hnd the minimum value of the following functional [8,9,10] 


See = 


4G 


d^y 


Vt 


Um 


1 + XL^R] +2XL 


'dM 


dyy hgK 


(3.3) 


where hij is the induced metric on the minimum surface, ho is the induced metric of 
the boundary of the minimum surface, R is the Ricci scalar of the surface and K is the 
extrinsic curvature of the extremal surface’s boundary. The same functional (3.3) can 
also be obtained using the derivation of gravitational entropy in [29,30,31,32], see also 
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[33,34], In the conformal case, where the gravitational dnal is the AdS 4 ^ space, it is easy to 
compnte the entanglement entropy of a disk of radins R. The minimnm snrface is described 
by r{z) = \/R‘^ — and the entanglement entropy is 



Let ns emphasise that, in general, there will be two kinds of constant terms in See, 
similar to those in (3.4). The hrst type of constant term, independent of the Ganss- 
Bonnet conpling, comes from the area of the minimnm snrface. However, this term is not 
topological as it receives corrections when the disc is deformed [35,36]. The constant term 
from the Ganss-Bonnet term, on the other hand, is proportional to the Enler characteristic 
and is topological by dehnition. As mentioned earlier, one of onr main goals is to hnd a 
model where only the constant term of a second kind is nonzero. 

The black hole entropy is no longer jnst the area of the horizon dne to the presence 
of the Ganss-Bonnet term. The general formnla for black hole entropy for the higher 
derivative gravity is the Wald entropy formnla [37] 

S = [ (fy\fh [ (fy\fh{l + XL^R). (3.5) 

N J horizon J horizon 

Here, hij and TZ are the indnced metric and Ricci scalar on the black hole horizon. The 
Lagrangian density C can be read off from the action (3.1). The binormal to the horizon 
e^i, is dehned as eab, where are two nnit normal vectors of the horizon 

and Cab is a nsnal Levi-Givita symbol. 

We are mostly interested in the constant term (the Enler characteristic) which is 
prodnced by integrating the Ricci scalar in (3.5) over the hypersnrface. To donble check 
onr prescription and to ensnre that no other constant terms are present, we can nse the 
observation of [38]. They showed that in a conformal theory, entanglement entropy of a 
disk eqnals the entropy of the hyperbolic black hole living in the dnal AdS space. To be 
precise, one can introdnce the coordinate transformation 

/z — pcosh('u) -I- \/— LX cosh(t/L), 

Lx^ / z — — L2sinh(t/L), 

(3.6) 

Lx /z — psinh('u) cos 6*, 

Lx^!z — psinh('u) sin 6* . 
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This coordinate transformation maps the metric of the empty AdS 


L2 


+ rjijdx^dx^) , 


(3.7) 


into the hyperbolic black hole with R x bonndary, with a horizon^ at p = L 


ds^ = 


dp^ 


—^ — 1 ) di + p^{du^ + sinh 'udd). 


p2/L2-l \L^ 

Using the Wald formnla (3.5), the entropy of the hyperbolic black hole (3.8) is 

7rp^(l — 2A) 


(3.8) 


A = 


2G 


du sinh('u) 


'u=0 


p=L 




(3.9) 


The valne of cosh •Umax = R/^ can be read off from the coordinate transformation (3.6) at 
p = L. Now we can see that the physical, cntoff-independent, terms in (3.4) and (3.9) are 
the same. This verihes the formnla for the entanglement entropy in (3.3). 


3.2. Gauss-Bonnet contributions to entanglement entropy and to the black hole entropy 
with genus g horizon 

The 2-dimensional Riemann snrface of genns g can be obtained by identifying the 
hyperbolic space by a hnite snbgronp of isometry [40] (see also [41] and refer¬ 
ence therein). In this case the Ganss-Bonnet term in (3.5) is proportional to the Enler 
characteristic, Xg, of the horizon dne to the Ganss-Bonnet theorem^ 

= —^(1 - d) , (3.10) 

where is the contribntion to the entropy dne to the Ganss-Bonnet term. Now, consider 
the Ganss-Bonnet term in the holographic entanglement entropy. The entangling snrface 
can be fonnd nsing the nsnal Ryn-Takayanagi prescription. In general, there are two types 

^ The appearance of the hyperbolic black hole here is similar to the way the Rindler space 
appears once the Rindler coordinates are used to describe the wedge of the Minkowski spacetime. 
See e.g. [39] for more detailed explanations. 

^ In the presence of boundaries this formula is modihed; in particular, the Euler characteristic 
of a disk is Xdisk = 1- 



of surfaces that extremize the area functional, the ones with the cylinder topology with 
r{z oo) = const and those with the disc topology with r{z — zq) = 0, for some hnite 
zq. The Euler characteristic is zero for the cylinder and unity for the disc. Let us denote 
the Gauss-Bonnet contribution in the entanglement entropy by . For the minimum 
surface with the disc topology, 5'^^ can be written as 

sS = ^ / Aviii? = (3,11) 

where we used Xdisk = 1- The expressions (3.10) and (3.11) are topological since they only 
depend on the topology of the horizon and minimum surface. The relation between (3.10) 
and (3.11) is the same as (1.3). 


4. A soft-wall holographic model with nonvanishing topological entanglement 
entropy 

In this section, we consider the geometry that is soft-wall model [42] in the IR and show 
that it satishes the criteria (i) and (ii) stated in the introduction. In the Appendix C, we 
also analyze the bottom-up model for fractional quantum Hall effect [43] in a certain range 
of parameters. Without the Gauss-Bonnet term in the action, the topological entanglement 
entropy is zero. However, since (i) and (ii) are satished, these models may have nonzero 
topological term when the Gauss-Bonnet term is included. 

The condition (i) will follow from the equation of motion, unaffected by Gauss-Bonnet 
term. We will see that there is no minimum surface with cylinder topology for the large 
values of the radius R. To ensure the condition (ii), one has to End the minimum surface 
and calculate the entanglement entropy to see that there is no 0{R^) term in the area of 
the minimum surface when R is large. 

We consider a class of geometries that are related to the AdS 4 by a warp-factor a{z): 

ds"^ = ^ {dz‘^ — dt^ -I- dr^ -\- r^d9^) . (4.1) 

The equation of motion can be obtained by minimizing the area of the surface described 
by r{z): 

A = I (4.2) 
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The equation describing the minimum surface is 




dz \z‘^y/l + r'(z)^ I 


(4.3) 


The metric is constructed to have a crossover scale between the UV and IR geometries set 
by the mass scale, //. In the UV region, nz 1, the warp factor is chosen to be 


auv{z -)■ 0) = 1, o!uv{^ ^ 0) = 0. 


(4.4) 


The second condition in (4.4) is chosen for technical convenience, so that the minimum 
surface of a disc radius R near the boundary, z/L <^l\s the same as in AdS^: 

.v2 


r[z) = ^ - 0(l/fl") 


(4.6) 


In the IR region, iiz ^ 1, the warp factor is chosen to be the so ft-wall warp factor of 


[42]. 


amiz) = e 




(4.6) 


The low energy spectrum of this theory can be found in existing works on soft-wall models. 
For = 1, the spectrum has a gap, set by the energy scale n, and continuous spectrum 
above the gap [44], For u > 1, the spectrum becomes gapped and discrete and for z/ < 1, 
the spectrum becomes gapless (see [45,46] and references therein). Note that, although 
the computations in the literature on soft-wall models are done in 5-dimensional gravity, 
the same conclusion can be reached in 4 dimensions, following the discussion in appendix 
A. The high energy spectrum of bound states in these models is given by 


2 2—2lv 

mt =n ’ 


(4.7) 


where n is the excitation number of the bound states. For z/ = 2, one obtains a behaviour 
similar to the Regge trajectory in QCD [47]. We will focus on the gapped system, z/ > 1, 
where the topological entanglement entropy is well dehned. 

As an explicit example, the warp factor that satishes the above condtions is 

a{z) = 1/cosh ^a/-I- 1 — 1 j (4.8) 

It is clear that the warp factor (4.8) reduces to the IR form (4.6) a,s iiz ^ 1 (with the 
substitution \/2L—)-L). In the UV regime, /uz 1, the warp factor above behaves as 

a{z)z:il-^{^zY'' + ... (4.9) 
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satisfying conditions in (4.4) for all the range of v where the spectrnm is gapped, 1 < u . 
The prohle of the snrface with the warp factor in (4.8) can be fonnd nnmerically and the 
relevant nnmerical resnlts are shown in Fignre 1. 

We will now show that the snrface with disc topology is the only possible solntion when 
the snrface probes the IR region, as long as z/ < 2. The eqnation of motion describing the 
minimnm snrface in the IR can be written as 





d [r(z)e r'(z) \ 

dz y a/1 + r'(z)2 J 


(4.10) 


We will now employ the method ontlined in the appendix C of [48] to rnle ont the solntion 
with cylinder topology. For the cylinder topology solntion to exist, there mnst be a solntion 
to (4.10) with the asymptotic solntion, r(z —)■ cxd) = cq. Hence, the cylinder solntion, at 
large z, mnst behave like 

r(z) = co + ciz^ + ... , (4.11) 


with m < 0. Note that (...) denotes the snb leading term in ^ cxd limit. Plngging in the 
ansatz (4.11) into the l.h.s. of (4.10) and extracting the leading term, one hnds that 


. 2^2 


z 2 


+ ... . 


(4.12) 


Similarly to the r.h.s. of (4.10), one hnds that 


d [r{z)r'{z)e \ 

dz y + r'{zY J 


= mcociz^ (—i'(juz)‘'— 3 + m) 

+ mc2z2m-4g-(/4^)- _ 3 + 2m) + ... . 


(4.13) 


We can see that, for the leading term on the left and right of (4.10) to match, we need the 
power of 2 in these terms to be identical, namely 


^-2 ^ ^i.+m-4 


z/ + m = 2 . 


(4.14) 


However, for z/ < 2, we can see that (4.14) cannot be satished. Therefore, an extremal 
snrface with the cylinder topology is not allowed for u < 2. Interestingly, this indicates 
that the phase transition for a disc region occnrs precisely at z/ = 2, the point of linear 
conhnement. On the other hand, for the slab region (see Appendix B), the phase transition 
occnrs at z/ = 1, where the spectrnm changes from gapped to gapless. 
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We now proceed to show that the constant term in the large radins expansion of 
the entanglement entropy vanishes. This can be done by showing that the area of the 
minimnm snrface does not contain the term, where R is the radins of the disc region 
at the bonndary. We will consider the case of z/ < 2 where only the snrface with the disc 
topology is a solntion. The tip of the snrface is located ai z — zq and iizq 3> 1 at large R. 
The mass gap, n, is set to nnity for the rest of this section (eqnivalently, all dimensionfnl 
qnantities are measnred in the nnits of //). Following the method ontlined in [49], we split 
the minimnm snrface into three parts: 

(I) Deep UV region, e < z < zi^^ , where e is the UV cntoff: This region contains part of 
the snrface that attached to the circle of radins R at the bonndary. The npper limit, 
zi^\ is the crossover scale where a{z) change the behavionr from auv{z) to am^z). 
As z 1, the minimnm snrface is described by (4.5) as demanded by the constrnction 
of a{z) in (4.4). 

(II) Intermediate region, zi^'^ < z < Zc^^ : In this region, the soft-wall warp factor, a(z), 

/o\ 

becomes ajji(z) in (4.6). The npper limit, Zc , is chosen snch that the area of the 
snrface in this region is not exponentially snppressed by the warp factor a/ji bnt will 
be snppressed when z > Zc^^. It is possible to hnd the prohle of the minimal snrface in 
the deep interior of this region i.e. when 2 : 3> 1 bnt z zq. To do this, we introdnce 
a new coordinate u = zjzQ and r{z) = {zq)'^p{u). The power of ( 2 : 0 )"" is chosen snch 
that p{u) is of order ( 2 : 0 )^. The eqnation of motion (4.10) in this new parametrisation 
becomes 


0 = 


(^0) 


2 —2n 


-lyu'" ^{zo)’'p'{u)- 


pin) 


P'{u) 

pin) 


— uu 




'(^0) 


[p'{u)Y + p" [u] 


(4.15) 

Collecting leading terms in large zq expansion, one hnds that the smallest valne of n 
that gives a nontrivial eqnation of motion is n = 1 — z//2. For this valne of n, the 
snrface in this region is described by 


P = y (^('^ - 2)di - ^2 '^), (4.16) 

where di is an integration constant. Expanding the solntion at small u and writing it 
in the original parametrisation, one hnds that 


r{z) 






(4,17) 


(21 

(III) Deep IR region, z > Zc where 2 can be of the same order as 20 : as 2 becomes very 
large, the area of the snrface in this region is exponentially snppressed by the warp 
factor aiR^z). 
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Fig. 1: (LEFT) Illustration of the region (I), (II) and (III). (RIGHT) Numerical 
value of radius R versus the position of the tip zo with a{z)~^ = cosh{\/z'^'^ + 1 — 1). 
The gradient of this plot is 1 — i^/2 for i/ = l(blue), 1.2 (orange) and 1.4 (green), 
respectively. 


From the solution r{z) in region (I) and (II), one can see that the minimum surface 
can be described by the following large R expansion, 

= (4.18) 

where we identify R ~ The relation between R and zq also agrees with the 

numerical results in hgure 1. Here {'f'i{z)} are functions interpolating between region 
(I) and (II). Note that this expansion breaks down when u approaches the critical value 
u = 2. This is expected from the previous analysis since the minimal surface might change 
its topology at this critical value of u. Plugging the expansion (4.17), into the area of the 
minimum surface (4.2), one hnds 


A 

A(u) 


^(I) + ^(II) + ^(III)5 
.^( 2 ) 

27rL^ R ' 




.( 1 ) 






a{z) 




+ 0{1/R^) 


5 


(4.19) 

where ^(i), ^(ii), ^(iii) correspond to the area of regions (I),(II) and (HI), respectively. In 

the region (I) the solution is approximately a cylinder of radius R in the large zq limit, 

and yields a typical UV divergence A cx L‘^R/e . The area of region (III) is exponentially 

suppressed by construction and can be neglected. To compute ^(ii), we hrst note that the 

hrst cross over scale is of the order 1/fi i.e. 4" ~ 1. Moreover, due to the fact that area 

(2) (2) 

from the region z > Zc is negligible, the upper limit, Zc can be lifted to inhnity without 


13 

















drastically changing the integral • Thns, hnite part of the area of the minimnm snrface 
can be written as 


^finite ~ ‘2,71 L 



1 


R 




n{z)^+0{l/R^) 


(4.20) 


and one can show that all integrals are hnite and independent of R. This indicates that 
there is no constant term, R^, in the area of the minimnm snrface. 

One can also check the validity of the above approximation scheme by compnting the 
entanglement nnmerically nsing the warp factor (4.8). We fonnd that the R^ term in the 
area (4.2) has a valne of order 10“^, for 1 < iz <2 and iiR ~ 20. This n um erical valne is 
negligibly small even thongh we are not at the limit /iR 3> 1 and is expected to decrease 
even fnrther as fiR increases. 

Upon inclnding the Ganss-Bonnet term, the entanglement entropy can be written as 

See = cxiR + + 0(1/R) , (4.21) 


where ai is some non-nniversal constant depending on the cntoff of the theory. Note that 
the bonndary term in (3.3) only gives a divergent term and therefore does not affect the R- 
independent piece. We can read off the topological entanglement entropy 7 by comparing 
(4.21) to (1.1), 


7 = 


ttLW 

G 


(4.22) 


5. Ground state degeneracy for the soft-wall model 

In this section, we make an attempt at compnting the gronnd state entropy for the 
model considered in the previons section when the horizon is the two-dimensional hyper¬ 
bolic space. As mentioned earlier, the snrface of genns > 1 is obtained by identifying 
this hyperbolic space by a finite snbgronp of the isometry. As in the empty AdS 4 , 
the area of the horizon is non-vanishing even at T = 0. However, this area is snppressed 
when the prodnct of the mass gap and the size, L, of T,g is large^. 

^ It would be interesting to do the same calculation for the fractional quantum hall model of 
[43]. Unfortunately, unlike the flat horizon case, we find no hyperbolic black hole solution when 
the dilaton has a scaling form </> ~ in the deep IR region for the allowed value of 7 and s 
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To find the black hole entropy with the horizon being one can proceed as the 

following. First, we replace the fiat spatial metric by the metric. The black hole 

metric has the following form 

= + ^ + . (5.1) 

where the spatial part of the bonndary, dTi\, is the line element on H^. 

[de^ + smh^ed(p'^) . (5.2) 


As mentioned earlier, the black hole entropy contains a contribntion from the Ganss- 
Bonnet term, even if the Ganss-Bonnet term is not dynamical. The entropy in this case 
can be compnted nsing (3.5). 


*5 = 


1 fL‘^a{zH), 


AG 


-'H 


vol(Sg) ) + -^Xh , 


(5.3) 


where the first term on the most right hand side is the area of the horizon and xh is the 
Enler characteristic of the genns g horizon i.e. xh = 2(1 — ( 7 ). 

In the following, we follow [50], where a semi-qnantitative method to constrnct the 
black hole solntions in the soft-wall models is proposed. In this setnp, the dilaton is 
assnmed to be non-dynamical and does not affect the metric in the string frame 

dz^ 


ds^ ■ 

string 


^2 


-f{z)dt^ + j^^+dEl 


ds^ = e 


(5.4) 


The dilaton ip{z) is chosen snch that it takes the form ip{z) — —{gzY{2, as in [47,50]. The 
emblackening factor is farther assnmed to be that of AdS-Schwarzschild 

f{z) = l-^ + MzY (5.5) 


This solntion is assnmed to be a solntion of a certain gravity model. To onr knowledge, 
snch model has not been fonnd^. 

® Given the metric of the form ds^ = + dx^ + dy^ + dz^), one can try to use the 

potential reconstruction method [51] to find a dilaton potential V{(f)) when A{z) = {gzY /2 + \ogz. 
The attempt to find F((/)) for black hole metric, ds^ = f {z)dt^ + dx^ + dy^ + dz^/f{z)), 

can be found in [52]. However, the dilaton potential for black hole phase is temperature dependent, 
so this route does not work for us. Another way to construct the soft-wall black hole is considered 
in [53] for a flat boundary and in the limit where the horizon is close to the boundary. In this 
method, the potential, V{(p,T), made out of dilaton, (f){z), and an additional scalar field, T{z), 
is reconstructed from non-black hole geometry [42]. To find a black hole solution, one has to 
solve for a nontrivial profile of f {z), A{z), 4’{z) and T{z). Unfortunately, we are unable to find an 
extremal hyperbolic black hole solution using this method. 
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The zero temperature solution can be found by tuning M into M^xt — (2/3^/^) 
and the horizon is located at zh = \/3-b. As a result, the warp factor at the horizon can 
be written as 


L'^aizn) 








(5.6) 

Substitute the expression in (5.6) back into (5.3), one hnds that the area of the horizon is 
exponentially suppressed when fxL 3> 1 . 

Hence, in this limit, the entropy of the soft-wall theory on Eg x contains only a 
Gauss-Bonnet term, which is topological 


S 


3 


27rAL2 

G 


(1-9). 


(5,7) 


Relating the expression in (5.7) to the topological entanglement entropy is (4.22), we hnd 
the relation (1.3). We emphasize that this result should be taken with a grain of salt, as 
a number of assumptions has been made to arrive at (5.7). It would be great to hnd an 
honest way of constructing holographic gapped geometries with hyperbolic horizons. 


6. Discussion 

We show that it is possible to obtain nonvanishing topological entanglement entropy, 

7 , in holography. The Gauss-Bonnet term plays a crucial role in our construction since 7 
is proportional to the Gauss-Bonnet coupling. The key property for the entangling surface 
to have a disk topology in the bulk is satished by the soft-wall models we consider. It is 
interesting that the soft and hard wall models of conhnement are clearly distinct from the 
point of view of our work. It would be interesting to identify a held-theoretic reason for 
this distinction. 

Let us emphasize that the dehnition of 7 involves computing the entanglement entropy 
of a disk in a theory on a plane. This is contrasted with a different measurement of the 
topological order: the degeneracy of the ground state in the same theory compactihed 
on a genus g surface. We observe that the relation between the topological entanglement 
entropy and the contribution to the ground state degeneracy from the Gauss-Bonnet term 
strongly resembles the same relation for the Ghern-Simons theory. It would be nice to 
understand this better. 

To describe the holographic dual of the soft-wall model on a higher genus surface one 
needs to consider the bulk action that leads to the soft-wall metric in the infrared and hnd 
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the solution with the asymptotic boundary being x R. The analogous procedure in 
the conformal case (pure AdS) leads to the asymptotically AdS black hole with hyperbolic 
horizon, so the appearance of the horizon would not be surprising. Unfortunately we did 
not succeed in constructing an honest soft-wall solution with a hyperbolic horizon due 
to technical difficulties. However, we present some arguments which indicate that the 
horizon area (and therefore contribution of the Einstein-Hilbert term to the ground state 
degeneracy) is exponentially suppressed as the mass gap becomes large. It would be nice 
to make these arguments more precise. 

To summarize, in this paper we showed how in certain holographic models with 
Einstein-Hilbert and Gauss-Bonnet terms, the held theoretic degrees of freedom dual to 
the former are frozen in the infrared due to conhning geometry, while the latter presumably 
give rise to a topological theory. To better understand the nature of this theory, more work 
is needed. In particular, considering holography on spaces with boundary can provide new 
insights. 
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Appendix A. Spectrum of gauge invariant mode in soft-wall model 

In this section, we extract gauge invariant combinations of the metric fluctuations 
in the soft-wall model and argue that the spectrum of the metric fluctuations in the 4- 
dimensional model is gapped as in the 5-dimensional one. 

The spectrum of the metric fluctuation of the 5-dimensional soft-wall model with gen¬ 
eral u was studied in [45]. Nevertheless, the metric fluctuations in 5 and 4 dimensions are 
slightly different. The 5-dimensional metric fluctuations, can be categorised 
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into three different channels by the remaining 0(2) symmetry of the bonndary field the¬ 
ory (see e.g. [54]). In this case, it is enongh to consider the scalar channel, consisting 
of only and solve scalar field eqnation of motion, where x* denote the directions 

along the bonndary. On the other hand, in 4-dimensional gravity, the finctnation is cate¬ 
gorised by the parity, y —)■ —y (see e.g. [55]). As a resnlt, the odd parity channel contains 
(hyt, hxy, hzy). The rest of the metric and scalar fields finctnations are in the even parity 
channel. The finctnations are conpled to components in their respective channel and, a 
priori, it is not clear that the analysis in 5-dimensional theory is still valid. 

We follow the approach similar to those in [55] by finding a gange invariant combina¬ 
tion of the finctnation in odd parity channel. The eqnations of motion of these finctnations 
are 


(a;2 _ P)hzy(z) - ikh' (z) + ujh'{z) = 0 , 




-a(z)^/^(h' (z) - ikhzy(z)) 

z 

-a(zfl‘^ {h'yt(z) + iujhzy(z)) 


-F u (uhxy(z) + khyt(z)) = 0, 
- k (khyt(z) -t- uhxy(z)) = 0. 


(A.l) 


where a{z) — e~O0 jg the warp factor of the soft-wall metric (4.1) in the IR. These three 
eqnations of motion are not independent. We can rearrange the last two eqnations into 
the eqnation of motion of the gange invariant combination by eliminating hzy(z). The 
resnlting eqnation is 




^ ! M/2 n ^ 
—a[z) ' y) (z) 


k 

+ -k'^)if{z) =0 , y:){z) = Ky{z) +-hyt(z) . (A.2) 

U) 


The eqnation of motion above is identical those of the KK modes in 5-dimensional soft-wall 
[46]. (see also [44] and [45] for discnssions in z/ = 1 and v = 2 case). Thns, we conclnde 
that the spectrnm of the metric finctnations in the 4-dimensional soft-wall is also gapped 
as in the 5-dimensional one. 


Appendix B. Entanglement entropy for slab geometry 

In this section, we stndy the entanglement entropy for the slab region in the field 
theory dnal to the model described by the metric (4.1). The case of soft-wall model with 
V = 2 was briefly mentioned in [56] and fnrther stndied in [57]. See also [58], for mnltiple 
slab regions in different gap phases. We fonnd that the phase transition of the entanglement 
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entropy occurs at z/ = 1, precisely the value when the spectrum changes from being gapped 
to gapless [45,46]. 

Let us introduce the notation here. The slab region is the region between x = ±£/2 
and has an inhnite length along y—direction. The induced metric on the minimum surface 
is written as 

^ ^2^^ ((1 + z\x)^)dx^ + dy'^) (B.l) 

To be consistent with the main text, the numerical result in this appendix is done with 
the warp factor a{z) — l/cosh(\/z 2 i^ _|_ \ where we set the energy gap scale = 1 

for simplicity. We will hrst look at the prohle of the minimum surface, determined by the 
area functional. The area of the surface in this section can be written as 

j-oo r^/2 

A = dy dx^/H{z)^/l + z'{xY ; H{z) = a{zY/z"^ , (B.2) 

J— oo J—il2 


As pointed out in previous studies, there are two possible conhgurations that satisfy 
the equation of motion derived from (B.2). The hrst solution is two inhnitely long parallel 
planes described by x{z) — ±^/2, referred to as disconnected surfaces. The other solution 
is the connected surface where z'{x) = 0 at some value of zq- To hud the prohle of z(x) 
describing the latter solution, we notice that area (B.2) does not explicitly depend of x 
and one can use the “conservation of energy” to obtain the following hrst order equation 


Inverting this equation of motion, one can obtain the relation between zq and the width £ 
of the slab as 


A/2 


r-Zo 


£^ 


dx = 


dz- 




(B.4) 


J-l /2 Jo \JH{z) — H{zq) 

The prohle of £ in (B.4) as a function of zq is shown in hgure 2. We found that for z/ > 1, 
the width £(zo) has a maximum value, £ — ^max- This indicates that there can only be the 
disconnected surfaces when £ > £max- The maximum value £ is not found for v < 1. 
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l(Zo) 



Fig. 2 : The width of the extremal surface l{zo) as a function of I/zq. The width 
i{zo) has a maximum value at large zq for n > 1. The maxima oi i{z) is not found 
for V < 1. This plot is obtained numerically from (B.4) 


To find the entanglement entropy, one needs to calculate the area (B.2) for both solu¬ 
tions and finds out which one is smaller. It turns out that the Gauss-Bonnet contribution 
for both solutions are zero and one can obtain the entanglement entropy by simply comput¬ 
ing the area. The area of both surfaces can be found by numerically evaluate the following 
integrals 


r-Zo 


= L^ 


dz 


H{z) 


^discon 


^H{z) - H{zo) 


= L^ 


dz^H{'. 


(B.5) 


Here e denotes the short distance cutoh while and represent the areas divided by 

length along y—direction of connected and disconnected surface respectively. The diherence 
between the areas of connected and disconnected surfaces, As = 5 “"^ — f = 1 

and 1.5 are shown in figure 3. We can see that, for u = 1, As approaches zero from 
below as we increase the width of the slab. Hence, for a slab region with a finite width, the 
connected surface remains a preferable solution. This is also true for the theory with z/ < 1. 
For u > 1, the minimum surface at small zq is the connected surface but undergoes the 
phase transition into disconnected surfaces at large zq, as depicted in figure 3 for = 1.5. 

To sum up, the entanglement entropy for a slab region with a large width ju£ ^ 1 is 
governed by disconnected surfaces for zv > 1 when the theory is gapped and by connected 
surface for z^ < 1 when the theory is gapless. As in [11], the critical point 1 / = 1,where the 
phase transition occurs does not ahect the entanglement entropy of a disk. This seems to 
indicate that the simple model of [56] actually needs more work. 
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- 1.0 


v=^.5 


'discon, 


'(i/=i) 


0.0 


0.1 


0.2 

1/Zo 


0.3 


0.4 


Fig. 3 : The difference of the areas, As, of the surfaces described by (B.5) in 
the unit of L^as a function of I/zq. We can se that for u = 1.5, there is a region 
where As > 0 indicated that the entanglement entropy is governed by disconnected 
surfaces when zo is large. For u = 1, the connected surface gives the smaller area 
for any zo- 

Appendix C. Bottom-up model of fractional quantum Hall system 

In this section, we apply onr procedure to a bottom-up model of quantum Hall system 
[43]. We show that this model has non-zero topological entanglement entropy at certain 
values of free parameters ( 7 , s) in the action of [43] when the Gauss-Bonnet term is present®. 

Let us briefly review the setup for this model. The action we consider here is the 
extension of Einstein-Maxwell-dilaton-axion theory, with four following components in the 
action 


(C.l) 


I = Ig + If + Iv + Igb , 


where Iqb is the usual Gauss-Bonnet term dehned in (3.1). Note again, that Iqb is does 


not affect the analysis in [43] and can be added to produce the topological entanglement 
entropy. The pieces Ig, If and ly can be written as 




(G.2) 


® In [43] the two free parameters are denoted as 7 and s. We denoted their 7 by 7 to a avoid 
confusion with the topological entanglement entropy 
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where V{4>,ri) is the SL(2, Z) invariant potential of the axion ri and dilaton cj). The 
solution we are interested in has both electric and magnetic helds however the electric 
held is completely screened out in the deep IR. In this region, the potential ri) is 
approximately V = —2Ae~^^^ and the extremal solution can be written as hyperscaling 
violating Lifshitz geometry 


ds‘^ = (7(7, s, h)p^ 


dV 


^2z 


+ 


dp^ + + dy'^ 


(C.3) 


p.. p- 

The overall constant, ( 7 ( 7 , s, h), depends on the free parameters in the action and the total 
magnetic held h. The dynamical exponent, 7, and hyperscaling violating exponent, 9, can 
be written in terms of 7 and s as 

(C,4) 

7 ^(l + s)(l —s) s — 1 

where p ^ 00 corresponds to the AdS boundary. We can dehne a new coordinate p = 

, X 2/6> 

( so that the boundary is at ^ 0 when 6 > 2. The induced metric on 

the entangling surface in this new coordinate is 


- — 
“'’ind 2 


1 + 


z'(r)^ 


dr^ _l_ 7*2^02 


n = 


29 

9^ 


(C.5) 


The allowed value of 9 (or equivalently 7 and s) can be found by imposing consistency 
conditions on the potential ti) and by demanding that the theory is gapped and has 
no naked singularity at hnite temperature. The results in [43] show that the allowed values 
of 7 and s are the small region around the line s = 1 ± 1 . 44(7 1) between 7 G ±[0.75,1]. 

Since the width in 7 direction is small, we approximate the allowed region of ( 7 , s) to be 
a straight line depicted in hgure 4. 



Fig. 4: (LEFT) The allowed value of s and 7 captured by the relation s = 
1± 1 . 44 ( 7 ± 1) where 7 G ±[0.75, 1] (RIGHT) The allowed value of n = 20/(6'- 2) 
from value for allowed value of 7 and s. Noted that n = 2 is the minimum value 
when s = 1 and 7 = ±1 
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The entanglement entropy calcnlated from the indnced metric of the form (C.5) has 
been explored in [48]. One of their resnlts is that there is only a solntion with a disc 
topology when n < 2. Thns, from condition (i), only theory with 7 = ±1 and s = 1 can 
have non-vanishing topological entanglement entropy. 

Moreover, when n = 2, the approximate solntion of the entangling snrface can be 
fonnd. Let ns parametrise the extremal snrface hy z — z{r), the eqnation of motion can 
be written as 


d [ rz'{r) A 2r{z{r)‘^ + 3z'{r)‘^) 

dr y 2 (r)^A/l + {z'(r)/z{r)) j z{r)^^/l + {z'{r)/z{r)y 


We are interested in the limit where the crossover from the deep IR geometry to the fnll 
metric happens ai z — Zc{Rc) and that Rc 3> 1, in the nnit of AdS radins. In the limit 
where r approaches zero, one hnds that the ansatz z{r) ~ solves the eqnation of 

motion (C. 6 ) at the leading order. After imposing the matching condition, z{Rc) = 1, to 
hx the constants Z and A, we hnd that the snrface in the region z > Zc is described by 


z{r) = exp 



(C.7) 


Now, we need to see whether there is a nonzero constant term from the area of the extremal 
snrface or not. We assnme that the entangling snrface extend very deep in the IR so that 
the hnite part of the area is determined by the IR part, similar to examples in section 4 
and in [49]. Also, in the large R limit, the crossover radins Rc is approximately eqnal to 
R, where R is the radins of the snrface at the AdS bonndary. In the large R zz R^ limit, 
we have 


Srt = 



(C. 8 ) 


This indicates that there is no constant term from this geometry. Hence, in this approxima¬ 
tion, this model passes the criteria (i) and (ii) mentioned in the introdnction. In principle, 
one shonld also extract the constant term in the UV completed metric, not jnst the IR 
part. However, the UV completed metric in [43] has to be obtained nnmerically which is 
beyond the scope of this work. 
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